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^ [ Abstract 

j^ ' Exact solutions of the perturbations equations in the magnetized ACDM scenario are 

presented. They apply during the dark ages and, more specifically, after the baryons are 
freed from the drag of the photons. The magnetized growth rate of matter perturbations 
v/-^ ■ is compared with the growth index obtained in the concordance paradigm and under the 

Tj- ■ of reionization to the stage of dark-energy dominance. The constraints derived from this 

f^ ■ analysis of are shown to be qualitatively complementary and quantitatively competitive with 

the bounds stemming from the analysis of the distortion patterns induced by the magnetized 
adiabatic mode on the temperature and polarization anisotropies of the Cosmic Microwave 
Background. 
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1 Introduction 

After the baryons are freed from the Compton drag of the photonqj (i.e. Zdrag = 1020.3 ±1.4) 
large-scale magnetic fields contribute to the evolution of the inhomogeneities in the spatial 
curvature and in the matter distribution even if observational constraints on their effects 
are lacking. After the drag epoch (i.e. z < -Zdrag) and before the Universe is reionized (i.e. 
z > z,:e = 10.5 ± 1.2) the Universe enters a very interesting epoch (often dubbed as dark 
age) when the first gravitationally bound system are formed. The purpose of the present 
paper is to explore the evolution of large-scale magnetic fields during the dark ages in the 
minimal framework of the magnetized ACDM scenario (niACDM in what follows) and in 
the simplest magnetohydrodynamical description which neglects the propagation of all the 
electromagnetic excitations of the plasma. 

Probably the best evidence of reionization comes from the analysis of the large-scale 
temperature and polarization fluctuations in the Cosmic Microwave background (CMB). 
Roughly speaking WMAP finds that about 8% of the CMB photons were scattered by 
free electrons in the intergalactic medium, but only 4 % could have been scattered by the 
intergalactic medium for redshifts z < 6. This means that the Universe must have been 
already reionized for z > 0(10). The direct estimates of Zj-e in terms of he WMAP 7yr data 
are consistent with the Sloan Digital Sky Survey (SDSS) jH [S] reporting the observations of 
quasars whose absorption spectra show, for z > 6, a substantial increase in the fraction of 
ionized hydrogen. In |6l [7] it has been argued, among other things, that the thickness of the 
last scattering surface and the optical depth at reionization are only mildly sensitive to the 
presence of an ambient magnetic field; the effect of reionization has then been parametrized, 
as in the conventional situation, with a second Gaussian peak of the visibility function 
centered around the reionization time. Conversely the inhomogeneities of the geometry and 
of the sources are affected by the presence of the ambient magnetic field whose effects on 
the baryons and on the geometry can never be completely switched off. 

Prior to photon decoupling the electrons, ions and photons are tightly coupled together; 
semi-analytic and numerical treatments of this problem are available in the presence of 
large-scale magnetic fields O [7] . Useful applications of the analytic treatments contemplate, 
for instance, explicit formulas for the distortions and for the shifts in the positions of the 
acoustic peaks as well as of the polarization peaks. For z < Zdrag the evolution of the in- 
homonogeneities in the mACDM scenario simplifies in many respects and the evolution of 
the large-scale inhomogeneities can be solved exactly. The solutions for the evolution of the 
metric inhomogeneities in the mACDM paradigm can also be extended to complementary 
cases provided the dark energy does not cluster as it can be conveniently assumed for differ- 
ent applications. In the literature there exist some exact solutions of the evolution equations 
of metric perturbations in pressureless media like the ones of Hwang [121 [13] or the ones 



^In what follows the values of the redshifts z will be estimated in terms of the standard ACDM paradigm 
analyzed in the light of the WMAP 7yr data alone [H El [3]. 



discussed in Refs. [HI [TSl [IS] • In the present paper after deriving the evolution equations of 
metric inhomogeneities for z < Zdrag, the exact solutions of the perturbation equations will 
be reported in the presence of a fully inhomogeneous magnetic field in the magnetohydrody- 
namical approximation stipulating, as customarily done, that the magnetic field, the Ohmic 
current and the electric field are all solenoidal (comoving) vectors. 

It is relevant to mention that various projects or radio arrays aim at a direct scrutiny of 
the dark ages like LOFAR (low- frequency array) |8], MWA (Murchison Wide-Field Array) 
[9], the enhanced VLA (Very Large Array) [10] and the SKA (square kilometer array) \TT\. 
One of the common strategies of these programs is the direct observational tomography of 
the (redshifted) 21 cm emission. The logic is to construct radio arrays capable of mapping 
the three-dimensional distribution of primeval hydrogen. Such a spatial distribution with its 
structures, its dumpiness, possibly its polarization bears the mark of the evolution through 
the dark ages and, most notably, of the first galaxies and stars producing ultraviolet emission 
and hence holes in the distribution of the protohydrogen. The role of cosmic magnetism in 
this epoch should be more carefully explored but the growth of matter inhomogeneities 
during the dark ages is certainly affected by the presence of large-scale magnetic fields, as 
general physical considerations suggest and as the solutions presented hereunder confirm. 
In section [2] we are going to present the evolution of the inhomogeneities in the mACDM 
model and for z < Zdrag- In section [3] the exact solutions of the system will be presented. 
In section H] the exact solutions are applied to the explicit calculation of the magnetized 
growth rate which is then compared to the conventional result of the concordance scenario 
with the purpose of deriving en explicit bound valid over the inhomogeneity scales relevant 
to structure formation. 

2 Evolution equations in real space 

The drag redshift corresponds, approximately, to the moment when the baryons are freed 
from the Compton drag of the photons and it can be expressed in terms of the present critical 
fraction of matter and of baryons ^7\ : 

1291 uj^Q n I ^ , ,921 
1 + 0.659 wg^^o' 
gi = 0.313u;Mr''[l + 0.607 u;°i^o'1, g2 = 0.238 c^S^^q'', (2-2) 

where cumo = ^o^mo and, more generally, Uxo = h^Qxa for a generic species x. For the 
present purposes, it is practical to decompose the synchronous gauge fluctuations directly in 
real space as 

VXg,j{x, t) = a^{T)l^d,djh{x, r) + G^didj^x, r) - U^.V^^x, r)] |, (2.3) 

where we have assumed a conformally flat background geometry characterized by the scale 
factor a{T) as it happens in the case of the ACDM paradigm. The density contrasts of dark 
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matter, baryons and dark energy will be denoted, respectively, by 6c, S^ and 64^- The Hamil- 
tonian constraint stemming from the (00) component of the perturbed Einstein equations in 
the gauge (Q ifl 
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E and B denote, respectively, the comoving electric and magnetic fields. The momentum 
constraint is instead given by 



drE = -n^ii + w^c)n^c0de + n 
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+ |a2F(f,r). (2.6) 



In the magnetohydrodynamical approximation adopted in the minimal mACDM model the 
displacement current is neglected since the frequencies involved by the present considerations 
are very small in comparison with the plasma frequency and the typical length-scales much 
larger than the Debye length. In the latter approximation not only the comoving magnetic 
field, the comoving electric field and the comoving Ohmic current are all solenoidal, i.e. 



V-5 = 0, V-E = 0, 



V-J = —V-{VxB) = 0. 
An 



(2.7) 



The solenoidal nature of the Ohmic current is consistent with the absence of the displacement 
current, as shown in the last relation of Eq. (12. 7p where J has been expressed, using the 
Maxwell equations for the comoving fields, as V x i?/(47r). When the ionization fraction 
Xe drops almost suddenly from 1 to about 10~^u;Mo/i^bo the concentration of the free charge 
carriers diminishes but the drop of the free charge carriers is faster than the drop of the 
temperature so that the plasma parameter decreases 



5'pla 



VoUoXe 



24:6^ 



/C(3) 



TT 



y/XeTfbo = 2.3 X 10" 



0.02258 



1/2 



(2.^ 



where Vd = 47rA£)/3 is the volume of the Debye sphere. Ad is the Debye length, r^bo is the 
ratio between the baryon and the photon concentrations and ({3) = 1.202. The smallness of 
5'piasma guarantees the validity of the plasma approximation; moreover, since the masses of 
the charge carriers are much larger than the kinetic temperature of the corresponding species, 
the conductivity a is proportional to the inverse of f^piasma- This means that the electric fields 
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will be suppressed as inverse powers of the conductivity or as positive powers of the plasma 
parameter. Thus the electric fields can be dropped from the dynamical equations for h and 
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where Wde denotes the barotropic index for dark energy. The evolution equation for the 
baryons can be written as 



a^ph 2 



(2.11) 



where pb is the baryonic mass density. The corresponding equations for the dark-matter 
species are: 

dr9c + n9^ = 0, drSc = ^drh - 9c. (2.12) 

The evolution equations for the dark energy are instead given by: 
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where Ade = V'^de and 6de = V^6'de; Cde denotes the sound speed of dark energy which is 
assigned independently of the barotropic index. 



3 Analytic solutions of system 

It is practical to start the integration of the system from the evolution equation of the baryon 
velocity reported in Eq. (12. lip . By using, as integration variable, a = a/a^e, the result of the 
integration o Eq. (12. lip can be expressed in terms of conventional hypergeometric functions 
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where a* and 9h{x, a*) denote, respectively, the initial integration time and the corresponding 
(inhomogeneous) value of the velocity field at a*; F[a, b; c; z] = 2-^1 [a, b; c; z] denotes the 



hypergeometric function fIE{ \W\ . The scale factor itself (appearing, for instance, in Eq. 
fl2.3p . is normalized to 1 at the present epoch (i.e. oq = 1) so that the following relations 
hold between Ode and the corresponding critical fractions of matter and dark energy: 

O-de = --r = [ 7=; ) ) ^deO = 1 ~ ^MO; ^MO = ^cO + ^bO- (3-2) 

^de + J- \ i 'MO / 

It is well known [20] that, in the synchronous gauge, the coordinate system is only partially 
fixed. The remaining gauge freedom must be removed from the initial conditions to avoid the 
potentially dangerous presence of spurious (i.e. gauge) modes. This problem is handled by 
removing the remaining gauge freedom or by expressing the results in terms of appropriate 
gauge- invariant variables |2T1 [22] (see also [23] ) . It is customary to select the dark matter 
rest frame and choose 9c{x,a) = to completely fix the coordinate system. In the latter 
case Eq. (12. 6p reads: 

dE ST-L r f coho 
da 2 a 
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Awmo/ 

where Q^c = 1 — ^m and the contribution of the radiation has been neglected since the 
drag redshift occurs when the background is already dominated by matter. If the dark 
energy is incompressible, its fluctuations can be neglected (i.e. 9de = and 6dc = 0) but the 
contribution on the background obviously persists since 

na) = ifo /^' / i^m, fiM(«) = ^^]^. (3.4) 



Note that the relations of Eq. (13. 4p can be used for a swift derivation of Eq. (13. ip . Dark 
energy is incompressible in the context of the mACDM scenario (where Wde — > — 1 and, 
according to Eq. (I2.14p . ^^e = exactly) but also in all the situations where Wde can be 
considered close to —1 (where, approximately, ^^e — 0). If the dark energy is compressible 
(i.e. ^do 7^ 0) then c^g must be different from Wde since i^de < —1/3. The bounds on c^g 
are currently rather loose and we shall assume, as customarily done [23], that < c^g < 1 
in the context of the so-called wCDM scenariqj. Since the barotropic index Wde and the 
sound speed Cde are assigned indipendently the total pressure fluctuation inherits a non- 
adiabatic contribution which is proportional to {wdc — c^e) ESj and this is the rationale of 
the appearance of such a term in Eq. (12.140 . For sake of simplicity we shall choose Cde = 
and, in this case, Eq. (13. 3p can be directly integrated and the result is: 

E{x,a) = E{x,a^) +M^{x,Zde,Wde)J^l{a,Wde) 

+ Q^{Zde,Wde)SB{x,a)J'2i(y,U;de) -V^{x,Zde,Wde)J^3i0i,Wde), (3.5) 



^By analyzing different data sets in tlie liglit of the tcCDM scenario, the error bars on Wdo either increase 
or they are restricted to an interval pinning down the ACDM value Wdc = —1- Values Wdo < — 1 will be 
excluded since, in these cases, the background may evolve towards a singularity in the future. This is only 
a practical choice since, in principle the present discussion can also be extended to the situation Wdc < — 1 
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where J\f({x, z^e, Wde), Q^l^de, u)dc) and V^{x, 2;de, u^de) are defined as: 
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The functions J-i(a,Wde) (with i = I, ...3) are either conventional or generahzed hypergeo- 
metric functions [T8l [T9l: 
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where pFg[ai, ..., ap]bi, ..., 6^; z] denotes the generahzed hypergeometric function [THl fT9] 
with: 
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Concerning the results of Eqs. (13. 5p . (13. 6p and (13. 7p few technical comments are on order. 
The integrals required for the actual solution can be performed by making extensive use 
of the transformation formulas of the hypergeometric functions [18J . A particularly useful 
formula is the one stipulating that 

F[a, b, ; c; -/^^^^^c] = (l + /J^^'^^V " F[c - a, c - b, ; c; -/S^'"'^^]. (3.9) 

The repeated use of Eq. (13.91) simplifies various expressions which can be reduced, in some 
cases, to the following indefinite integral 
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which is a consequence of the formula giving the nth derivative of a generalized hypergeo- 
metric function [27]. The solution of Eqs. (13. 5p and (I3.6p -f l3^ solves also Eq. (I2.10p whose 
explicit form in terms of the notation used in this section is: 
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By knowing explicitly the solution in terms of S(x, a) it is immediate to compute the curva- 
ture perturbations on comoving orthogonal hypersurfaces whose explicit expression in terms 
of ^ can be written as: 



which also implies that 
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Since the curvature perturbations are invariant for infinitesimal gauge transformations they 
will have the same value in any gauge. In the CDM rest frame, adopted throughout to get 
rid of spurious gauge modes, 6c{x, a) can be immediately computed 
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Similarly the solution for 6b{x, a) turns out to be 
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The evolution equations of the the total matter density contrast (5m (x, a) 

5„ 
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can be derived by combining the previous equations. Defining the differential operator 
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the evolution of Sra{x, a) is given by: 
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Note that in Eqs. (13.171) and (13.181) the dark energy component has been assumed to be 
smooth and incompressible. As noted before this assumption holds exactly when Wdc = — 1 
(as in the mACDM case) and approximately when Wde — —1- The solution of Eqs. (I3.17p - 
(I3.18P can be written as 
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denotes the Wronskian of the two independent solutions of the homogeneous equation whose 
exphcit expressions are 
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Finally the solution for h{x, a) is given by 
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To normalize properly the solutions obtained for a > a,= it is useful, even if not mandatory, to 
compute analytically the solutions for a < a* by assuming that the dark energy background 
is negligible. In this case the matter-radiation system can also be solved in the presence of 
large-scale magnetic fields and we shall be interested in computing the obtained solutions in 
the limit a^, < a <^ Oeq- Defining, in analogy with Eq. (I3.16p . the differential operator C2{<y) 
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the evolution of 5^{x, a) in the range a^. < a <^ Oeq can be written as: 
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The solution of the previous equation reads 

(5m (x, a) = Vi{x)J^s{a) + V2{x)T9{a) + 
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where a^ = fli/cteq denotes the initial reference value of the scale factor while G^'^\a, (3, ) and 
W{f3) are defined as 

J^s{/3)J'9{a) - Ml3)Ma) . 1 



G^'\a,(3, 



WiP,w^e 



(3.28) 



The two linearly independent solutions J-8(a) and J^g{a) are given, this time, by the following 
expressions: 
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Following the remarks already proposed around Eq. (13. lip it is useful to reiterate that 
the gauge dependent description pursued here can be complemented by appropriate gauge- 
invariant treatments. For instance the variables 
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are gauge- invariant and become, in the uniform curvature gauge [121 [13], the density contrasts 
of the single species. In equally correct terms we could also argue that the (^ of the species 
X measures the curvature perturbations on the hypersurfaces where the energy density of 
X is uniform. From the perspective of the synchronous gauge the latter statements can be 
appreciated by noticing that ( is proportional to C, which is related to TZ (see Eq. (13. lip ) 
corresponding, in turn, to the curvature perturbation in the comoving orthogonal gauge. 

The results derived in the present section can be used for different purposes. In what 
follows the attention will be focused on the calculation of the growth factor of matter in- 
homogeneities and on its comparison with the results obtained in the more conventional 
situation where magnetic fields are assumed to be absent. 



4 Applications and discussions 

The integral appearing in Eq. (I3.27P can be evaluated explicitly and the obtained results 
expanded in the limit aj = a-Jacq > 1 and a = a/acq ^ 1 which corresponds to the 
physically interesting situation where the initial conditions of the solution (I3.19P are set well 
after matter-radiation equality and anyway for redshifts smaller than ^drag- In the latter 
limit the result can be expressed as 

.^L^ + I),(x)(l + -a)+0(-)+0(- 
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where, we remind, a = a/aeq^. While in the regime a,, < a <^ Oeq the total matter density 
contrast is well estimated by Eq. (14. ip , for a > a* the solution is instead given by Eq. fl3.19p . 
The initial conditions for the calculation of the growth rate can therefore be set by requiring 
that 

(4.2) 



d 
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da 



The two conditions reported in Eq. (14. 2 p fix the constants appearing in Eq. (I3.19p . The 
growth rate [261 [2H] (see also [29]) can be computed from Eq. (I3.19P and it is given by 
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where, according to Eq. (I3.19p . 5m(x, a) can be written as 
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The two constants Ci(wde, ^de, -z*) and €2(104^, Zde, z^) are a lengthy combination of hyper- 
geometric functions depending on a* = l/{z^ + 1) and a^e = l/(-2de + !)• The quantity 
6{x) denotes the spatial profile of the fluctuation of the total density contrast which can be 
estimated in terms of the matter power spectrum while Sb(x, a) denotes the contribution of 
the magnetic sources. The power spectrum of 6{x) can be written as 



{S{x)S{x + r}) = fdlnkPsik) 



sinfcr 



kr 
4 . / k\"--^, 2 



Vs{k,y,^) = ^^^[y) ln^(^/^eq), (4.5) 

holding for wavenumbers larger than fceq = 0.00974'^q'}5oq\qMpc^^. In Eq. (14. 5 p k^ = 
0.002 Mpc"^ and An = (2.43 ± 0.11) x 10"^ in the case of the WMAP 7yr data alone 
[H [21 13] and for the case of the ACDM scenario. 

Typical wavenumbers k ^ k^.^ crossed inside the Hubble volume before matter-radiation 
equality. Conversely the very large length-scales (relevant for the region of the Sachs- 
Wolfe plateau and for the integrated Sachs- Wolfe effect) fall into the complementary regime 
k ^ fceq- The wavenumbers touched by the present considerations range from fcmin = 
0.01/ioMpc~^ to, approximately, fcmax = 0-ShoMpc~^. The range kmin < ^ < ^max in- 
cludes also the scale at which the spectrum becomes nonlinear, i.e. k^i > 0.09 /loMpc"^. 
The illustration of the analytical results will be given in terms of the following fiducial set 
of parameters determined on the basis of the WMAP 7yr data alone [H |2l |3]: 

{Qw, ficc^deo, ho, n„ ere) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088), (4.6) 

where e^e denotes the optical depth at reionization and Us is the spectral index of curvature 
perturbations. The power spectra contributing to T,b{x, a) can be solely expressed in terms 
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of ^b{x, a) and o"B(a;, a) by recalling that 
V-(JxS) 



Aa^p^ 



WVb-V^^b, fiB(x,«) 



3 n a^py ' 



9,9,ni^' = -p.VVb, (4.7) 



where Hg denotes the anisotropic stress of the magnetic fields. In the range /cmin < ^ < 
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Figure 1: The growth rate is illustrated for two different values of the barotropic index, i.e. 
Wde = —1 (full line) and Wdc = —0.6 (dashed line). In the plot at the left i?L = InG and 
k = 0.03 Mpc"^ In the plot at the right B^ = 10 nG and k = 0.01 Mpc"\ In both plots the 
magnetic spectral index has been chosen as ub = 1-5. 

fcmax; fhc tcrms containing the Ohmic current in the evolution equations (e.g. the term 
containing Sb{x, a) in 7^ (x, a)) dominates against the one containing just flBix, r) because 
of the presence of two supplementary spatial gradients. The same comment can be made for 
S(x, a). In Fig. [T]the growth rate is reported for two illustrative values of the barotropic 
index w^e = —1 (full line) and w^e = —0.6 (dashed line); in both plots of Fig. [H fc = 
0.03 Mpc~^. In Fourier space the power spectra of ^b and o"b will be denoted by Vn{k) and 
Va{k) and their explicit expression can be written as 



Vn{k)=SB{-^ 



J. \ 2(nB-l)+an 



/ ]^ \ 2(nB-l)+«n 

V„{k)=rB£B[j- 



{U 



an and a^j are the corresponding running of the spectral indices which are set to zero in the 
minimal magnetized ACDM scenario. The relative amplitude of the two power spectra at 
the magnetic pivot scale /cl (conventionally chosen to be 1 Mpc~ ) is controlled by tb. Since 
the magnetic fields are, themselves, stochastically distributed, the ensemble average of their 
Fourier modes obeys 

-> Itt'^ -> / k \ "B— 1 

{B,{k)B,ip)) = —P,^(k)PBik)5^'\k + ff), VBik) = AB[-^) , (4.9) 

where Pij{k) = {k'^Sij — kikj)/k'^ and ^b the spectral amplitude of the magnetic field. The 
spectral amphtude of the magnetic energy density £b depends upon the ultraviolet cut- 
off associated with diffusion damping (i.e. /cd), upon the infrared cut-off associate with 
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the (comoving) angular diameter distance at last scattering (i.e. Dp^{zi)) and also on the 
magnetic spectral index n-Q, i.e. [6l [7] 






nB + 29 

20(7 -ne: 
nB + 29 

20(7 -ne: 



^ (5-2nB)(83nB-473) /_^\ i-" 
2(7-nB)(nB + 29) Ua/ 
(riB - l)(87nB - 501) / k \ 5-2n,B 
^ (nB + 29)(7-nB) VA^ 



riB < 1, 
nB > 1, 



(4.10) 
(4.11) 



where, recalling that Da{zi) = 2dA{zi) / (Hoy/^Mo) , 

k^iz,) 2240rfA(^i) 



kA{zi) = 1/Da{zi), 



^a(^i) s/V^^TT 



/rm 



^1 \ , ,0.24 -0.11 



(4.12) 



In Eq. (I4.12P tri is the ratio of radiation to matter at zi which is the redshift of last- 
scattering and which can be determined analytically for typical values of the parameters 
close to the best-fit determined on the basis of the WMAP 7yr data and in the light of the 
ACDM paradigm: 



'^Rl 



91 



PR[Zi} 



4.15 X 10- 



Pm[Zi) 
0.0783(tJb) 



UJm 
-0.238 



[l + 39.5(wb)°-^^=']' 



■(- 

U03 
92 



zi = 1048[l + fiio^''][l+giujf^ 



0.560 



i + 2i.i{uj^y-^^' 



(4.13) 



where /i and /2 are given, respectively, by /i = 1.24 x 10^^, /2 = 0.738. When n-Q = 1 
the spectra for the energy density and for the Lorentz force are scale-invariant once the 
logarithmic divergence of the two-point function is appropriately subtracted. The amplitude 
^B can be traded for the magnetic field regularized over a spatial domain k^^, i.e. in the 
case riB > 1 



Sv 



4(7 -^b) 



Al 



3(nB-l)(5-2nB)(87rp, 



Av 



(27r 



in.R-l 



r[(nB-l)/2] 



Bl 



(4.14) 



and analogously in the ub < 1 case but with ^b = [(1 — ^b) /2](kA/ kt)^^^"'^^ Bf^. 

In the case of the minimal mACDM the relation of (nB, ^b) to (nB, -Bl) follows from Eqs. 
(l4.9p -f H?T^ . The magnetic energy density can be naturally referred to the energy density of 



the photon background so that Qb l 
of curvature perturbations: 



BI/{8ttp^) can be measured in units of the amplitude 



n 



BL 



A 



39.568 



n 



Bl 
nG 



T. 



70 



2.725 K 



A 



n 



2.41 X 10-9 



-1 



(4.15) 



The parameter space of the magnetized tyCDM models and of the magnetized ACDM sce- 
nario have been investigated, respectively, in |6] and [7]. In a frequentist approach, the 
boundaries of the confidence regions obtained in P [7] represent exclusion plots at 68.3 % 
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and 95.4% confidence level. When moving from the magnetized ACDM scenario to the 
magnetized wCDM model we have that the the parameters maximizing the likelihood get 
shifted to slightly larger valueCl 



('^B, -Sl)acdm 
(^B, -Bl)acdm 



;i.598, 3.156nG) -^ (ub, 5l)»cdm = (1.883, 4.982 nG), (4.16) 
;i.616, 3.218nG) -^ (nB,5L)»cDM = (1.913, 5.163 nG). (4.17) 



Even if the addition of a fluctuating dark energy background pins down systematically larger 
values of the magnetic field parameters, the results of [HI [7] will be used here just for a 
consistent illustration of the results. In Fig. [2] the growth rate is illustrated as a function of 
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Figure 2: In the plot at the left the growth rate is illustrated for a = 1 (i.e. a = a^e) 
for i?L = InG (full line) and 5l = 10 nG (dashed line). In the right plot a = 0.1 (i.e. 
a = 0.1 Ode) and the values of Bl are, respectively, -Bl = 10 nG (full line) and -Bl = 30 nG 
(dashed line). In both plots ne = 1.5. and k = 0.03 Mpc^^. 

Wde and for two different choices of i?L, i-e. InG (full line) and 10 nG (dashed line). The 
difference between the plot at the left and at the right is the value of the redshift: while in the 
plot at the left a = Ode (i-e. z = z^e) in the plot at the right a = 0.1 a^e and z > z^e- In Fig. 
|3] the growth rate is illustrated as a function of Sl for different choices of the wavenumber 
and of the barotropic index. The full and dashed lines correspond, in each plot of Fig. |3l 
to two different values of the magnetic spectral index, i.e. ne = 1.5 (full line) and tib = 1.2 
dashed line. 

In the absence of large-scale magnetic fields, the growth rate can be parametrized as 
fi1i{a) where 7 is the growth index which can be explicitly estimated as [SD] (see also [2B] ) 



liWde 



3wdc - 3 3 (1 - w;do)(l - 3wdc/2) 



6wde — 5 125 



;i - 6w7de/5)2 



e + 0{e' 



(4.18) 



''The difference between Eqs. (|4.16p and (|4.17p is tliat tfie parameters of Eq. (|4.16p are obtained from 
the analysis of the temperature autocorrelations while the parameters of Eq. (|4.17p are obtained by adding 
the data points of the cross-correlations between temperature and E-mode polarization [B] . 
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Figure 3: In both plots the growth rate is illustrated for a = Oeq as a function of the 
regularized magnetic field intensity -Bl- In the plot at the left Wde = —1 (and k = 0.03 Mpc^^) 
while in the plot at the right w^e = —0.6 (and k = 0.3Mpc~^). 

where e = 1 — Qm- The requirement that the magnetized growth rate of Eq. fl4.3p does not 
exceed the standard fit for the growth rate implies an upper bound on Sl which is illustrated 
in Fig. m The bound on Bj^ gets more stringent as k increases and for large redshift. In 
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Figure 4: The bound on Bl illustrated for n-Q = 1.6 (long dashed line), n-Q = 1.5 (full 
line) and n-Q = 1.4 (short dashed line). In the left plot k = 0.03 Mpc~^. In the right plot 
k = 0.3Mpc~^. In both plots a = Oeq 

Fig. Hla = Ode- The allowed values of -Bl for different values of w^c stay below each of the 
curves reported in the plots of Fig. HJ Recalling the definition of the comoving magnetic 
Jeans length [29] 



A 



n 



BJ 



G a^Ph 



1.90 X 10-2 



l^bO 



0.02258 



-1 



^ I "p^- 



Bl 



ina-^ Pb ' 



(4.19) 



the results of Fig. H] imply Abj < (9(Mpc). The considerations reported in the present 
analysis suggest that the study of large-scale magnetism during the dark ages provide further 
constraints on the evolution of large-scale magnetic fields. The quantitative features of the 
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obtained constraints have been shown to be both complementary and competitive with the 
bounds stemming from a direct analysis of potential distortions induced on the temperature 
and polarization anisotropies for a range of length scales shorter than the one relevant for 
the physics of the CMB [32]. 

If the origin of the large-scale magnetic fields is primordial (as opposed to astrophysical) 
it is plausible to expect the presence of magnetic fields in the primeval plasma also before 
the decoupling of radiation from matter. If the evolution of the large-scale magnetic fields 
follows the tenets of magnetohydrodynamics then the magnetic flux will be conserved across 
the last scattering and potential effects of the magnetic fields during the dark ages could help 
to decide on their origin and implications. In this perspective the SKA pT] program could 
provide decisive informations both from the foreseen full sky surveys of Faraday rotation and 
also from more detailed pictures of structure formation and reionization over a wide range 
of redshifts during the dark age. 
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